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A COMPARISON OF FAR IR AND HI AS REDDENING PREDICTORS AT HIGH GALACTIC LATITUDE 

J. E. G. Peek^2 

ABSTRACT 

Both the Galactic 21-cm line flux from neutral hydrogen (Hi) in interstellar medium and the far- 
infrared (FIR) emission from Galactic dust grains have been used to estimate the strength of Galactic 
reddening of distant sources. In this work we use a collection of uniform color distant galaxies as 
color standards to determine whether the Hi method or the FIR method is superior. We find that 
the two methods both produce reasonably accurate maps, but that both show significant errors as 
compared to the typical color of the background galaxies. We find that a mixture of the FIR and 
Hi maps in roughly a 2-to-l ratio is clearly superior to either map alone. We recommend that future 
reddening maps should use both sets of data, and that well-constructed FIR and Hi maps should both 
be vigorously pursued. 

Subject headings: astronomical databases: atlases, ISM: dust, extinction, galaxies: photometry 



1. INTRODUCTION 

The interstellar medium (ISM) is a very sooty material 
by the standards of Earth's atmosphere, suffused with 
obscuring astrophysical dust. Through emission, absorp- 
tion, and scattering, this dust acts to modify the light 
we receive from distant objects across almost all observ- 
able wavelengths, most notoriously acting to artificially 
redden objects in the optical. One way we combat this 
distortion of our observations is to construct maps of the 
surface density of dust on the sky, such that observa- 
tions can be corrected back to an unreddened state. As 
we push forward into an era of precision cosmology, it is 
especially important to have high accuracy maps at high 
Galactic latitude, where the bulk of cosmological investi- 
gations take place (e. g. Huterer, Cunha, & Fang, 2013, 
[Menard et al. 2010). 

In the past two main methods have been used to con- 
struct such reddening maps. The first is the Hi method. 
The 21-cm hyperfine transition of the hydrogen atom 
traces neutral gas throughout the Galactic ISM. At 
high latitudes the Hi ISM is largely optically th in 



^Hi < 1; with some exceptions e. g. Peek et al. 



2011b) 



such that the integrated 21-cm profile is a good proxy 
for the column density. Since the high Galactic latitude 
ISM has a largely uniform dust-to-gas ratio, and neutral 
gas dominates the column density, the 21-cm intensity 
can act as a proxy for dust column. This method was 
largely codified by |Burstein fc Heile sI (1978), and maps 
were presented in Burstein ^ Heiles (1982; BH) . A more 
modern data base of Hi columns is presented in [Kalberla 



et al. (2005; LAB), which has the advantages of being 
more coherently reduced, all sky, and having careful cor- 
rections for stray radiation, which can be a major source 
of error in Hi column density maps. 

The Hi method acted as the standard for high Galactic 
latitude extinctio n until the far infrar ed (FIR) method 
was developed by |Schlegel etaTI ( p98| SFD). In the FIR 
method the assumption is made that there is a single 
population of dust, both in terms of grain size distribu- 
tion and composition, and that therefore there is a di- 
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rect relationship between the FIR thermal emission from 
grains and their contribution to extinction and redden- 
ing. This relationship is subject to the variations in the 
temperature of the grains. To build an extinction map 
SFD made an all-sky map from the IRAS 100 jam data, 
as well as a dust temperature map from DIRBE. Com- 
bining these two, they presented an all- sky extinction 
map which has been the standard since its publication. 
The SFD map superseded the BH map due largely to its 
higher resolution, fidelity, completeness, and ease of use, 
rather than some known superiority of the FIR method 
over the Hi method. 

Both the Hi and FIR methods have inherent flaws that 
stem from their assumptions. In reality, the Hi column 
is not entirely optically thin, and thus the integrated 
21-cm flux is not a perfect proxy for neutral column den- 
sity. Furthermore, other phases of the ISM carry ob- 
scuring dust, including molecular and ionized phases. 
We also know that the dust-to-gas ratio is not a fixed 
quantity, even with these assumptions taken into account 



(Burstein & Heiles 
known issues with t 



19781). There are also a number of 
leTTR method. SFD note that vari- 



ation of grain temperature along the line of sight can 
generate significant errors. Grain size variation can also 
effect the results since FIR dust emission is biased toward 
larger grains, but extinction is biased toward smaller 
grains. Indeed, some authors have found dust emissivity 
is ~ 3 times higher in molec ular gas than in atomic gas 
due to grai n agglomeration (Ade et al. 2011 Stepnik 
|et al.| |2003| . There are also other sources of FIK emis 
sion, especially unresolved, distant galaxies, which can 
contamina te these maps (Yahata et al.^ 2007( Kashiwagi 



et al.||2012 i). 

THe modern era of digitized optical surveys has al- 
lowed us to test the accuracy of these maps in some 
detail, specifically, with the extremely consistent photo- 
metric calibration of the Sloan Digita l Sky Survey (SDSS; 
York et al.||2000||Padmanabhan et al.||2'Q08p . There has 



been a burst of activity using ''standard crayon" meth- 
ods; selections of large sets of objects corrected to have 
identical intrinsic colors such that variations from these 
colors must be due to effects from the intervening dust. 



Schlafiy et al. (2010) used photometric observations of 
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large groups of stars, [Peek fc Graves] ( |2Q1Q[ PGIO) used 
spectroscopic obs ervations of non- star forming galaxies, 
Jones et a l. (2011*) used spectroscopic ob servations of cool 
stars, and Schlafly & Finkbeiner (2011| SFll) used spec- 
troscopic observations of a broader range of stars. All of 
these explorations found that SFD is largely accurate in 
the optical, but that overall errors and spatial variability 
of the reddening are detectable. 

In this work we attempt to resolve whether the FIR 
method, the Hi method, or some combination make the 
most consistent reddening maps. To do this we employ 
the standard crayons developed in PGIO and compare 
maps of the residual colors of these galaxies to the FIR 
maps of SFD to the Hi column densities derived from 
LAB. To put the Hi and FIR methods on equal footing 
we work at a resolution lower than that of either survey, 
which also serves to aggregate enough PGIO galaxies to 
give statistically significant color variations from the ex- 
pectations of SFD. In ^we briefly summarize the meth- 
ods used to produce the low-resolution maps of Galaxy 
colors. Hi , and FIR. In ^we show the relative accuracy 
of the two methods, and a combination of the methods. 
We discuss these results and conclude in Q 

2. DATA AND METHODS 

Here we aim to generate three maps to compare the 
Hi method predictions and FIR method predictions to 
the observed color variation of galaxy colors on the sky. 
We would prefer to have Hi maps at a resolution compa- 
rable t o the SFD FIR ma ps, but the only such large-area 
maps ([Peek et al. 2011a) are neither complete nor cor- 



rected for stray radiation, and thus are at present unsuit 
able for such an analysis. Furthermore, to get to a high 
enough signal-to-noise to detect significant color varia- 
tions from predictions we must average the PGIO galax- 
ies over larger areas than even the LAB survey beam. We 
note that the highe r resolution GASS survey (McClure- 
Griffiths et al.|2009D would be a useful map for this anal- 
ysis, were it in the same hemisphere as the SDSS data, 
and we expect the EBHIS maps will also b e useful once 
comp lete and corrected for stray radiation (Winkel et al. 
2010). Given this state of affairs we employ the LAB sur 



vey column densities as our data set from which to gen- 
erate our Hi method maps, and SFD for the FIR method 
maps. We n ote that while a 14% c orrection to SFD was 
suggested by 'Schl afly fc Finkbeiner | ( |2011|) , it is not clear 
that this correction applies to the very low extinctions 
weare investigating at high latitude, and [Peek fc Schim>| 
novich (2013j PS 13) find no evidence for such a variation 
in either the optical or UV. We therefore use the SFD 
reddening predictions with no correction. Throughout 
this text we will refer to the E {B — V) reddening pre- 
dictions from SFD as Esfd and the predictions from 
LAB as Elab- 

The 151,637 PGIO galaxies were selected from the 
SDSS DR7 main galaxy sample data release to have no 
detected Oil or Ha, and thus be almost entirely non- 
starforming, passively evolving galaxies. The galaxies 
were K-corrected, corrected for an absolute magnitude- 
color relationship as a function of redshift, as well as 
corrected for very weak dependence of galaxy color on 
galaxy environmental density. The final variation in ^ — r 
galaxy color translates to a typical standard deviation in 
E{B -V) of 0.024 magnitudes. The details of this anal- 



ysis are available in PGIO. As in PS 13, we employ the 
galaxy colors as uncorrected by SFD. 

We follow the same map-making methods as described 
in PS 13. Essentially, we put the galaxies into a zenith 
equal area projection and grid them into 16 square de- 
gree pixels. To evaluate the E {B — V) in these pixels we 
use the median color of the galaxies within each pixel, as 
the color distribution has outliers that can corrupt the 
average color. We refer to this value as Epqiq- The 
error in each pixel of the PGIO map is evaluated using 
a standard boostrap method. Each pixel has a typical 
noise of 2 millimagnitudes in measured E {B — V). We 
restrict our investigation to areas of sky where we have 
errors in this map below 4 millimagnitudes. To provide a 
reliable comparison, we evaluate the SFD and LAB maps 
at the same positions as the galaxies within a given pixel 
and take the median value within the pixel. While in 
practice this makes little difference on our final results, 
it nicely preserves the weighting within each pixel that 
comes from the precise galaxy positions within the sur- 
vey. This is especially relevant towards the edges of the 
survey. 

3. RESULTS 

The resulting maps of EsfDi ^hi^lab^ and Epcio are 
shown in Figure [T] Also shown are the residuals when 
taking the difference between pairs of maps. To do this 
we perform a first-order polynomial error-weighted least- 
squares fit solving 



and 



EpGio — «o + clsfdEsfd 



EpGlO = ^0 + ^LAbNhFLAB' 



(1) 

(2) 



We find ao = -1.36 x 10"^, a^FD = 1-10 and 6o = 
-3.78 X 10"^,6lab = 1.44 x 10"^^. As this analysis is 
performed over a somewhat restricted area of sky we do 
not specifically recommend these values for future work, 
but rather report the m to show broad consist e ncy with 
previous work (e. [Burste in fc Heiles||1978| (|Schlafiy| 
fc Finkbeiner 2011)).^ As expected, the three Esfd^ 
J^Hi.LABi and EpGio maps look rather similar. We 
note 'the that the residual maps (PGIO-SFD, PGIO-LAB) 
have many pixels that are significantly different than the 
expectation, in a statisitical sense. What is perhaps most 
interesting is that the PGIO-SFD and PGIO-LAB resid- 
ual maps have inconsistent deviations from zero. This is 
to say, whatever processes make the Hi method fail are 
not the same processes that make the FIR method fail. 

To further quantify this discrepancy we analyze the 
map data in the PGIO-SFD versus PGIO-LAB space, 
as in Figure |2j We find the residual maps are cor- 
related with amplitude of ac = 0.0023 and the anti- 
correlated with amplitude da = 0.0017. The preponder- 
ance of the correlation stems from the noise in the PGIO 
data itself (which contributes to both maps) with mean 
c^PGio = 0.0020 (see We find that the amplitude of 
the correlation in the residual maps that stem from the 
HI and FIR maps themselves is only acr = 0.0010, less 
than 60% of the anticorrelation. We note that some fair 
fraction of acr may in fact be due to errors in the SDSS 
data or the PGIO sample which are not captured by our 
bootstrap error analysis. Imperfections in the SDSS pho- 
tometric calibrations, for instance, would contribute to 




Figure 1. Maps of reddening towards the northern Galactic cap. The top row shows the SFD, LAB, and PGIO maps described in ^ 
We do not show the smah fraction of sky covered toward the southern Galactic cap. The bottom row shows the difference between various 
maps. To take the difference we simply take the residual of a first-order polynomial error-weighted least-squares fit of the PGIO map to the 
SFD map and the PGIO data to the LAB map. The difference of these two fits is shown at the far right. We note that the color residuals 
from the Hi and LAB maps are both highly statisitically significant (with typical erorrs of 0.002) and are structurally different, implying 
that their errors stem from different effects in the ISM. 
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Figure 2. A scatter plot of the PGIO-SFD versus PGIO-LAB 
maps, with filled circles representing each pixel. In red is the 1-cr 
ellispe that shows the correlation (towards y + x) and anticorre- 
lation (towards y — x) between the two maps. The amplitudes of 
these dispersions is shown at the lower left of the figure. At the 
bottom center of the figure is the expected correlation between the 
maps that stems simply from the known errors in PGIO, crpGio- At 
bottom right is the residual correlation between the maps, which 
we note is significantly smaller than the anticorrelation. 

this error (see e. g. PGIO figure 7). Even without ap- 
pealing to sources of error beyond the statisitcal error 
in the PGIO map, the anticorrelation between the two 
residual maps is larger than the correlation, indicating 



that the FIR and HI maps are complementary. 

With this in mind we perform a similar error- weighted 
least-squares fit to the following equation: 

EpGlO = Co + CFIRasFoEsFD + CHlbLAsNni^LAB (3) 

using the previously determined values of asFD and 
bhAB^ and solving for cq^cfir^ and chi- Additionally, 
we run a simple bootstrap analysis to probe the error 
space, resampling the map pixels with replacement. The 
results are shown in Figure [Sj We find best fit values of 
cfir = 0.66^ Chi = 0.37. This to say a healthy mixture 
of the two methods, with a weighting toward the FIR 
method, produces a significantly superior map to using 
either map alone. Errors in the PGIO-SFD and PGIO- 
LAB maps are 0.0042 magnitudes in both cases, as ass- 
esed using the standard deviation of the data clipped at 
3.5 sigma. The equivalent error in the best fit mixture 
map is 0.0038 magnitudes. 

4. DISCUSSION & CONCLUSION 

Our analysis shows that while both the FIR method 
(in this case SFD) and the Hi method (in this case LAB) 
produce maps that are largely error- free, each map does 
have detectable variation from the observed reddening 
(PGIO). Furthermore, we show that these variations are 
not well-correlated with one another - whatever mecha- 
nism produces discrepancies in the Hi reddening map is 
not the same mechanism that produces discrepancies in 
the FIR reddening map. Quantitatively, an admixture of 
66% FIR and 37% HI outperforms either map alone. 

This has significant implications for the construction 
of future reddening maps. Specifically, it shows that 
the FIR method and Hi method have complementary 
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errors in each map, which we expect to remain despite 
improvements in resolution. The second is that each one 
of these methods may be individuahy improved, thus ren- 
dering the other less useful. The FIR methods could be 
improved with higher fidelity and resolution temperature 
maps, as well as corrections for extragalactic sources of 
FIR light. The Hi maps could be improved with some ad 
hoc correction for Hi opacity, and perhaps the addition 
of molecular tracers like CO or OH, or even corrections 
for ionized gas. While these kinds of advancements are 
likely to improve the quality of each method, their basic 
weaknesses as described in ^T]are unlikely to disappear, 
and so the combination of the two methods will still be 
superior to either one alone. 

JEGP would like to thank Doug Finkbeiner, Eddie 
Schlafly, David Schiminovich, Karin Sandstrom, and 
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Figure 3. The best fit to Equation |3] is shown as a red cross. 
Gray points are the results of 1000 Monte Carlo bootstraps, to 
give a sense of the errors. Clearly, a 66% FIR + 37% HI map is 
superior to either a 100% FIR or 100% Hi map. 
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